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We propose a simple way to measure the three-dimensional rotational diffusion of micrometric 
wires, using two-dimensional video microscopy. The out-of-plane Brownian motion of the wires in a 
viscous fluid is deduced from their projection on the focal plane of an optical microscope objective. 
An angular variable reflecting the out-of-plane motion, and satisfying a Langevin equation, is 
computed from the apparent wire length and its projected angular displacement. The rotational 
diffusion coefhcient of wires between 1 — 100 fim is extracted, as well as the diameter distribution. 
Translational and rotational diffusion were found to be in good agreement. This is a promising 
way to characterize soft visco-elastic materials, and probe the dimension of anisotropic objects. 

1 Introduction 

Recording the three-dimensional (3D) motion of anisotropically shaped probes is a challenging is- 
sue. Although the theory of rotational Brownian motion has been established for a long time now—, 
the direct visualization and quantification of the Brownian motion of a micrometric anisotropic 
probe with a microscope is recent. It is due to the difficulty of quantifying the three-dimensional 
motion of the probe with two-dimensional (2D) optical techniques. It was first solved by studying 
the in-plane rotational motion of anisotropic probes ^'3.. Recently, highly specialized optical tech- 
niques have opened new opportunities. Rotational diffusion was studied using light streak tracking 
of thin microdisks^i^, depolarized dynamic light scattering and epifiuorescence microscopy of opti- 
cally anisotropic spherical colloidal probes~, scanning confocal microscopy of colloidal rods with 
three-dimensional resolution^, or reconstruction of the wire position from its hologram observed 
on the focal plane of a microscope*. The rotation along the long axis was also investigated by 
analyzing the fluorescence images of rodlike tetramers^. Following the 3D rotational diffusion of 
an optical probe thus remains costly in equipment, as well as in computational power. 

From a practical point of view, micrometric probes can be used to determine the relation 
between stress and deformation in materials reducing significantly the sample volume, which may 
be crucial in biological samples^''. The technique, called microrheology, is a powerful tool to 
probe the rheological properties of complex fluids and biological materials at the micrometric 
scale. It can be achieved, either by recording the thermal fluctuations of probes immersed in the 
material, or by active manipulation of the probes^. While microrheology based on translational 
diffusion has been extensively investigated, the rotational diffusion of anisotropic objects remains 
poorly explored. However, it may be of great interest to investigate the length-scale dependent 
rheological properties of heterogeneous structured materials, such as complex fluids or biological 
tissues. In the case of anisotropic probes such as wires, the large aspect ratio of the probe allows 
for a detectable Brownian motion, over larger length scales, typically between 1 — 100 /im, than 
for spherical probes. 

In this letter, we propose a simple way to measure the 3D rotational diffusion of micrometric 
wires, using 2D video microscopy. The 3D rotational Brownian motion of the wires immersed in 
a viscous fluid is extracted from their 2D projection on the focal plane of a microscope objective. 
An angular variable reflecting the out-of-plane motion of the wires and satisfying a Langevin 
equation, was computed from the apparent wire length and its projected angular displacement. 
The rotational diffusion coefficient was found to vary over more than 5 decades, for wires of length 
between 1 — 100 /im and anisotropy ratios in the range 2 — 2000. The resolution of the technique was 



1 



Figure 1: Spherical coordinate system {r,6,ip) for a wire diffusing in a three-dimensional space. 
The unit vector Ur is the wire orientational vector. 



quantified, by analyzing the wires trajectories. From these measurements, we were able to extract 
the distribution of the wires diameter, in good agreement with electron microscopy measurements. 
Rotational and translational diffusion measurements were compared, giving good agreement^ii^. 
This provides a new step towards the reliable use of rotational diffusion to characterize complex 
materials with an optical microscope, and probe the dimension of anisotropic objects. 

2 Langevin equations for rotational diffusion 

Let us consider a rigid wire diffusing in a stationary viscous fluid. Its rotational Brownian motion 
can be modeled by a Langevin equation, describing the fluctuations of the wire orientational unit 
vector Ur (Fig.[T]). In absence of an external torque, and neglecting inertia, the rotational equation 
of motion writes : 



in which (^hlj and are two effects from the fluid, respectively the viscous drag and the random 
Langevin torque. Since the wire is axisymmetric, the matrix of friction coefficients (^r is diagonal 
in the frame {ur, Ue, u^). The eigenvalue along Ur describes the friction opposing the self- rotation 
of the wire around its main axis. The eigenvalues along iie and Uip are equal, and describe the 
friction opposing the rotation of the wire main axis. 

The projection of the rotational Langevin equation ([T]) on the plane perpendicular to the wire 
{ue, Uip), leads to : 



random torque, in the plane {ue, u^p), and Cr is the friction coefficient perpendicular to the wire 
axis. The projection of equation ([T]) along the wire axis will not be considered here^. 

The random torque can be written as T^^ = — Ti (t) ue + T2 [t) Uip, where Ti{t) and 
T2{t) are two Gaussian white- noise thermal driving torque, satisfying : 



= -Ci?w + . 
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with T the bath temperature, and ( ) refers to a time-averaged quantity. 
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For a finite cylinder - length L and diameter d - the perpendicular friction coefficient can be 
written in the form : 

where g{L/d) is a dimensionless function, which takes into account the finite-size effects of the 
wire. This function was analytically calculated for ellipsoids—, and numerically approximated in 
the case of cylinders^ii^""— . 

The rotation vector ijJ-^ can be expressed as a;-*- = Oii^ — sinO (pug, and the projection of 
Eq. [2]in the plane (ug,Uip) thus gives : 

Cr^ sine if = Ti{t), CR^0 = T2{t) 

The angular variables ip{t), defined such as -ip — sinO ip, and d{t) both obey a one-dimensional 
Langevin equation, which respectively leads to : 

(AV'2(i)) = 2^t = 2Dnt (4) 

{A0^t))=2DRt (5) 

where Dr is the rotational diffusion coefficient. Computing Eqs. (jH) and the diffusion coeffi- 
cient simply writes : 

Dr = '^9iL/d) (6) 

In the general case of an out-of-plane rotational diffusion, Eqs. ^ and ([5]) show that deter- 
mining the variables o^' ^(^) '^iU lead to the rotational diffusion coefficient Dr. From the 
2D video recordings, both (p{t) and smd{t) were extracted. The variable Aip{t) — ijj{t) — ijj{0) — 
J^dt' sin e{t') ifit') was then computed, leading to the determination of the mean-squared angular 
displacement {Atjj'^ {t)) , and therefore to the rotational diffusion coefficient Dr. 



3 Material and Methods 

The wire formation results from the electrostatic co-assembly between oppositely charged iron 
oxide nanoparticles and polymers^. The wires are purified and suspended in DI water. Figured] 
shows an optical microscopy image of the wire suspension after synthesis, and the corresponding 
length distribution in inset. In the present study, two batches of wires of length 15 fj,m (Fig. [2]) and 
30 fim are investigated—. Due to a rather broad polydispersity in length, wires of length between 
1 and 100 /Ltm were obtained. The distribution of the wires diameter was determined from electron 
microscopy, with a median diameter 0.4 /xm, leading to anisotropy ratios between 2 — 2000. 

The aqueous wire suspension was then mixed with pure glycerol giving aqueous solutions of 
glycerol, also referred as wire suspensions. Aqueous solutions of glycerol with two different volume 
fractions, 50% and 60%, were prepared. The wire suspension was then introduced in an observation 
chamber (3 mm x 3 mm x 250 /im) between a microscope slide and a coverslip, sealed with araldite 
glue to avoid evaporation and contamination of the sample. 

An inverted Leica DM IRB microscope with a xlOO oil immersion objective (NA=1.3, free 
working distance : 130 /im), coupled to a camera (EoSens Mikrotron) were used to record the 2D 
projection of the wires thermal fluctuations on the focal plane objective. The wires concentration 
was chosen diluted enough to prevent collisions and hydrodynamic coupling. They were always 
tracked far enough from the walls of the observation chamber. The microscope objective tem- 
perature was controlled within 0.1° C, using a Bioptechs heating ring coupled to a home-made 
cooling device. The sample temperature was controlled through the oil immersion in contact. 
Sedimentation of the wires was negligible on the recording time scales. 

The camera was typically recording 10 images per second during 200 s (2000 images). The 3D 
Brownian motion of the wires was extracted from their 2D projection on the (x, y) plane (Fig. [1]). 
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Figure 2: Optical microscopy images (X 40) of a wire suspension. The distribution of the lengthes 
after synthesis can be fitted by a log- normal distribution, as shown in inset (83 wires). Due to a 
rather broad polydispersity, of the order of 0.5 (defined as the standard deviation of Ini), wires 
of length between 1 to 100 fim were obtained. The distribution of the diameters was determined 
from electron microscopy, with a median diameter 0.4 /im, leading to anisotropy ratios between 
2 - 2000. 



The angle ip{t) and the projected length L{t) were measured from the images, using a home-made 
tracking algorithm, which is implemented as an ImageJ pluginl-^.. Since the objective depth of 
focus (1 — 2 /im) is shorter than the length of the wires, their out-of-plane image is distorted. The 
distorted image of a wire was seen as a cluster of connected pixels. The algorithm output is the 
length L{t) of the cluster (projected length of the wire) and the angle (p{t) of the cluster, as defined 
in figure [TJ The algorithm mainly consists in three steps. First, a user-defined threshold is applied 
to the image. Then, the wire projection, seen as a set of connected pixels (cluster) is tracked at 
time t, in the vicinity of the position at time t — 1. Finally, the cluster orientation and its length 
are computed, respectively giving the angle ip(t), and the projected length L(t). The position 
of the cluster center was also computed. Wires with a high out-of-plane angle, corresponding to 
angle 9 smaller than 50 degrees, could not be considered. Since the tracked wires are chosen to 
lie in the focal plane at the beginning of the recording, the length of the wire L was taken as the 
maximum measured length L{t) within the recording time, leading to sin 6 — L(t)/ max {L{t)). 

The quantity Aiplt) was computed by using the discrete equation Aip{t, to) = X]t''=*o~'^* 

where Stp{t') = (p{t' + 6t) — (f{t'), and 6t is the time lapse between two images. A time average 

then enables us to compute (Ai/)^(i)) — ((Ai/)^(t, to)))to- 

The uncertainty on (Ai/)^(t)) was evaluated including the uncertainties on Lp, the projected 
length, and the total length. It also takes into account statistical accuracy. The uncertainty on 
the apparent length of an out-of-plane wire was determined by varying the z-position of the focal 
plane while recording the fixed wire. The corresponding computed lengths for different z-positions 
of the focal plane give an uncertainty of the apparent length as the wire moves along the axial 
direction. It was estimated to be 8% for wires longer than 4 /im and 15% below. 

4 Rotational diffusion coefficient and diameters distribu- 
tion 

The mean-squared angular displacement (MSAD) {Atp'^{t)) is shown in Fig. [3] as a function of 
the lag time t, for wires of length between 3 /tm and 100 /tm. The MSAD {AijP'it)) was found to 
increase linearly with time, as expected in a viscous fluid. From these curves, a rotational diffusion 
coefficient Da defined such as (Ai/)^(i)) = 2D at, could then be extracted. Figure |4] shows the 
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Figure 3: Mean-squared angular displacement (MS AD) of the wires as a function of the lag 
time, for wires of length between 3 and 100 /im. The experiments were performed in an aqueous 
solution of glycerol (50% in volume). The MSAD increases linearly with time, as expected in a 
purely viscous fluid. The slope of each curve, corresponding to the diffusion coefficient, decreases 
with the length of the wire. Similar data were obtained in the 60% in volume aqueous solutions 
of glycerol (not shown). The error bars include the uncertainties on L, tp and statistical accuracy. 
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Figure 4: Rotational diffusion coefficient of the wires as a function of the wires length, normalised 
by viscosity and temperature. At first order, the diffusion coefficient scales as L~^, as expected 
from Eq. The dimensionless parameter grot{L/d) — §§^Dii is shown in inset as a function of 
the wires length. The gray line corresponds to the analytical expression established by Broersma, 
assuming a median diameter of 400 nm. The large distribution of the data around the gray line 
reflects the distribution of the wires diameter. 
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Figure 5: Measurement error trot as a function of the length of the wires. The error in orientation 
is around 4° for a 10 /im long wire. The largest error was obtained for short wires, of 2 to 3 
micrometers long, where an error of orientation of roughly 8° was found. For the longest wires 
above 20 /zm, the error was found to be less than 1°. 



rescaled quantity -^^^Dji as a function of the measured length of the wire L. At leading order, 
the diffusion coefficient decreases as L^^, over 5 decades. A correction grot{L/d) for the finite-size 



effects of the wire is expected, as described in Eq. The rescaled diffusion coefficient 



was then multiplied by , leading to the experimental determination of the dimensionless function 
Qrot (L/d) (Fig.Hinset). 

The resolution of the technique was quantified from the wires trajectories using the relation : 



(7) 



including the measurement error e^ot, and corrected for the camera exposure time a^^. The 
measurement error erot was found to decrease with the length of the wire, as shown in figure [5] It 
typically corresponds to an error in orientation of 4° for a 10 long wire. The largest error was 
obtained for short wires, of 2 to 3 micrometers long, where an error of orientation of roughly 8° 
was found. For the longest wires above 20 /im, the error was found to be less than 1°. 

An analytical expression for the finite-size effects function grot{L/d) was established by Broersma 
in 1960 ^^1^° . It is expected to be valid for p = L/d > 4.6, and writes : 



= In (p)- 0.446 
16 63 



0.2 

ln(2p) 
64 



\n{2py ln{2py \n{2py 



(8) 



Fig. ID-inset shows the Broesma relation ([8|), using a median diameter d = 400 nm. The large 
distribution of the data around this adjustment reflects the distribution of the diameter of the 
wires. 

By numerically inverting the relation ([8|), we could extract the distribution of the wires di- 
ameter. The corresponding values oi p — L/d range between 2 and 2000. A small fraction of 
the wires (less than 5%) exhibit a value of p which is smaller than the validity range established 
by Broersma {p > 4.6). More recently, another analytical expression, valid for 2 < p < 20, was 
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Figure 6: Distribution of the wires diameter, extracted from the wires Brownian motion measure- 
ments. The distributions are fitted by log-normal distributions, leading to equal median diameters 
0.4 /im. The polydispersities are respectively crdmicro = 0.53 and adsEM — 0.38. The distributions, 
obtained with differents methods are in good agreement. In inset, scanning electron microscopy 
of a wire. 



proposed by Tirado et al.^*. This expression was used for the few wires falling out of the validity 
range of Broersma's theory, and for which Tirado's expression is valid. 

Figure IHl shows the distribution of the diameters obtained from the diffusion measurements, 
also compared to the one obtained from Scanning Electron Microscopy (SEM) measurements. 
The distributions have been fitted by log-normal distributions, leading to equal median diameters 
0.4 /zm. The polydispersities (standard deviations of Ind) are respectively crdmicro = 0.53 and 
CdSEM = 0.38. The distributions are in good agreement, the largest discrepencies come for the 
smallest diameters, where the precision of our method is the lowest, because of the high non- 
linearity of the function grot{L/ d) as a function of d. 

We now compare rotational and translational fluctuations. The translational fluctuations of the 
wires were measured in the wire's frame of reference. In our experiments, only the component of 
the center-of-mass translation along u^p could be measured. Since the diffusion in the z-direction 
could not be evaluated, the translational diffusion parallel to the wire, along m^, and the one 
along iig were not accessible. Figure [7] shows the center-of-mass mean-squared displacement 
(Ar^(t)) = {{{r{t + t') - r{t')) ■ u^{t')f)t' for wires, between 3 and 100 ^m. It was found to 
increase linearly with the lag time, as expected in a viscous fluid. The data of translational 
diffusion were fitted according to the relation—: 

(Ar^(t)) = 2Z?i„„,,(t - a/3) + 2eL„.i (9) 

including the measurement error etransi and corrected for the camera's exposure time a. The 
translational diffusion coefficient writes : 

^ kgT _|_ 

D transl = /^^jj^StransliP) (10) 

with the finite-size effect function gi^^nsiip) = ^^(p) + 0.839 + 0.185/p + 0.233/p2ii. The measure- 
ment error etransi is shown in Figure [7] (inset) for wires of length between 3 and 100 /xm. It was 
found to be less than etransi — 0.3 pixel in most cases, which means that the center of the wire 
was tracked with a precision better than 40 nm in the (x, y) plane. 
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L(/ini) 


Db (radVs) 






In(ptransi) 


4.4 


0.016 ±0.004 


39 ±4 


3.3 ± 1.1 


3.4 ±0.6 


9.8 


(1.3 ±0.2)10-^ 


17.9 ±0.7 


3.2 ±0.7 


3.6 ±0.3 


35 


(4.9 ±0.2)10-^ 


6.7±0.2 


5.2 ±0.7 


5.2 ±0.2 


95 


(2.8±0.1)10~« 


3.1 ±0.2 


5.8 ±0.7 


6.8 ±0.5 



Table 1: Comparison between rotational and translational measurements for wires of various 
lengths. The values of ln(p) obtained with rotational and translational measurements, respectively 
In(prot) and ln{ptransi) , are in good agreement. 



Fitting Eqs (O and ([9]) for a given wire gives the values Djf and D^^^^^^i of the diffusion 
coefficients, respectively obtained with rotational and translational measurements. Given the ex- 
perimental conditions, this respectively yields values of ln(p), In(prot) and lTi{ptransi), according 
to Eqs. ([ni and ([TU]). Table 1 presents a comparison between rotational and translational mea- 
surements for wires of length between 4 and 95 /xm. The values of ln(p) obtained with rotational 
and translational measurements were found to be in good agreement. 

5 Conclusion 

In this letter, we propose a simple way to follow the 3D rotational Brownian motion of micrometric 
wires in a viscous fluid, from the 2D projection of the wires on the focal plane of a microscope. 
The rotational diffusion coefficient of the wires between 1 — 100 /im was computed as a function 
of the wires length, ranging over 5 decades. The resolution of the technique was quantified, by 
analyzing the wires trajectories. Our diffusion measurements allow us to extract the distribution 
of the wires diameter, in good agreement with SEM measurements. Rotational and translational 
diffusion measurements were compared and found to be in good agreement This technique 
provides a simple way to measure the out-of-plane rotational diffusion of a wire in a viscous fluid, 
opening new opportunities in microrheology to characterize more complex fluids, and probe the 
dimension of anisotropic objects. 
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